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Research Article 

A Laboratory Study on the Molecular 
Basis of Primary Congenital Glaucoma 


Abstract 

Purpose: To detect pathogenic mutations in cytochrome P450 familyl subfamily B polypeptidel 
(CYP1B1) gene in nineteen sporadic Primary congenital glaucoma (PCG) cases and to identify patients 
lacking CYP1B1 mutations. 

Methods: CYP1B1 exon 2 and the coding part of exon 3 of 1 5 participants were amplified by 
Polymerase chain reaction and amplicons were sequenced by Sanger sequencing. Sequencing data was 
analyzed to identify the gene mutations or Single Nucleotide Polymorphisms SNPs. 

Results: Four previously reported PCG-associated CYP1B1 mutations (c. 1159G>A; p.E387K, 
c.230T>C; p.L77P, c.l 103G>A; p.R368H and c.l 568G>A; p.R523K) were found in four patients out of the 
15 fully 'sequenced' patients. Also, ten previously reported Single Nucleotide Polymorphisms and two 
novel noncoding variants were identified. 

Conclusion: The relatively low percentage of PCG patients having CYP1B1 mutations (4/15=26.6%) 
demonstrates that other known and unknown genes may contribute to PCG pathogenesis. Lack of CYP1B1 
gene mutations in some patients stresses the need to identify other responsible candidates. 


Introduction 

Primary congenital glaucoma (PCG) is a classical form of 
infant buphthalmos and a predominant type of congenital 
glaucoma [1,2]. Although, buphthalmos has been reported 
from the time of Hippocrates (460-377 BC), it was associated 
with high intraocular pressure (IOP) and vision loss in children 
in the middle of the eighteenth century [1,3,4]. In PCG, isolated 
maldevelopment of ocular drainage structures increases 
resistance to aqueous outflow [3,5,6]. This leads to increased 
IOP, optic nerve atrophy and blindness if prompt and proper 
diagnosis and management is not provided [4,7,8]. 

Infants with PCG are affected differently. However, 
presentation of epiphora, blepharospasm, photophobia and 
swollen or cloudy cornea due to high IOP, is the norm [7]. 

Familial, like sporadic, PCG is an autosomal recessive 
disease [9]. It is highly prevalent in the consanguineous and 
inbred Gypsy subpopulation of Slovakia. This was reported to 
be due to genetic drift or ‘founder effect’ [8,10]. In non Gypsy 
populations, PCG may be inherited in a multifactorial fashion 

[11]. 

Although the genetic component of PCG is debatable, 4 
chromosomal loci has been mapped including, GLAUCOMA 


3, PRIMARY CONGENITAL, A (GLC3A on 2p22.2 locus; OMIM 
231300), GLC3B (ip36.2-p36.l; OMIM 600975), GLC3C (l4q24 
.2; OMIM 601771) and GLC3D on I4q24.3 locus (OMIM 613086). 

Recently, PCG phenotype was associated with mutations in 
cytochrome P450 family lsubfamily B polypeptide 1 ( CYP1B1 ) 
gene on the 2p2i locus, gene symbol GLC3A (OMIM #231300) 
and Latent transforming growth factor beta-binding protein 2 
( LTBP2 ) gene (GLC3D; OMIM 602091) [8]. The genes associated 
with GLC3B and GLC3C are not yet known [12]. 

CYP1B1 gene analysis is the focus of this study. CYPiBigene 
(Ensembl transcript ESNT00000260630) is a three-exon gene 
in which only exons 2 and part of 3 are translated [9]. It encodes 
a 543-amino-acid ‘drug metabolizing’ protein or enzyme 
involved in early anterior chamber angle development [13]. 
Mutations in CYP1B1 have also been implicated in other anterior 
segment dysgenesis (ASD) syndromes [14]. The molecular 
role of CYP1B1 gene in the pathophysiology of PCG is not fully 
understood [8]. This study identifies unrelated PCG subjects 
that have or lack CYP1B1 mutations by screening the coding 
exons of the CYPiBigene using direct sequencing procedures. 
Patients having no disease-causing CYP1B1 alleles were noted 
for exome sequencing to find out other genetic factors involved 
in PCG pathogenesis [9]. 
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Materials and Methods 

Participants 

Informed consent was obtained from the parents or 
responsible guardians of the subjects. The probands and their 
immediate families were questioned on their past medical 
history and examined at the hospital. Reviewing the eye charts 
from their respective eye physicians helped confirm their visual 
function. The inclusion criteria were intraocular pressure 
>2immHg in the first years of life, the presence of cornea 
oedema, scar, or Haab striae or an ocular history consistent 
with infantile glaucoma. The presence of enlarged globe in their 
teens or adulthood suggests a history of congenital glaucoma 
in the first few years of life. Other clinical signs were examined 
using portable slit lamp and direct ophthalmoscope. 

Congenital glaucoma CYP1B1 mutations (Table 1) 

The study did not report the individual visual acuity, visual 
field and intraocular pressure of the patients. 

DNA quantification and quality assessment 

In order to determine the concentration and quality 
of extracted DNA (patient samples), the Nanodrop 8000 
(Thermo Scientific) was used to obtain the absorbance of the 
DNA samples at 26onm; A 26o =i was noted as the equivalent 
to 50ng/ml for double stranded DNA. The Nanodrop machine 
was ‘blanked’ initially by adding of nuclease-free watenpl 
(dH 20 ) onto the tip at the base of the spectrophotometer. The 
concentration and 260/280 values were noted (not shown). 


Table 1 : Congenital Glaucoma CYP1 B1 Mutations (JH, TC, HZ, MF, BY are patient's 
name code, lab no. Labouratory number. DOB; Date of Birth, Homo; Homozygous, 
Het; heterozygous) 

Patient 

Mutation 

Homo/Het 

JH DOB 20/8/90 MG55312 
(lab no. 11005694 & 067246) 

P.E387K (rs55989760) 

Homo 

TC DOB 25/12/10 Ml 1 /003053 
(lab no. 11000358) 

p.R368H (rs28936414) 
p.L77P(CM000136) 

Het 

Het 

HZ DOB 1 6/9/08 MG69428 
(lab no. 10004412) 

p.R368H (rs28936414) 

Homo 

MF DOB 7/11/08 E08/1 6439 
(lab no. 089365) 

p.E387K(rs55989760) 

Homo 

BY DOB 9/3/04 SP04/3021 
(lab no. 041104) 

p.R523T (CM066028) 

Het (missing 2 nd 
mutation) 


A 260/280 absorbance ratio of > 1.8 was considered a better 
amount of purity for DNA samples. Samples with concentration 
of more than loong/pl were diluted properly with dH 20 to 
make loong/pl concentration. 

The DNA samples (in loong/ul concentrations) were spun 
down for 3-5 secs using Minispin at a speed of 2000rpm to 
prevent DNA from being shattered. 

Designing primers 

In order to design primers of the CYP1B1 gene (DNA) 
target, the following databases were used; ensembl; a 
massive database/genome browser, primer3plus and reverse 
complement. 

Human CYP1B1 (ENST00000260630) gene sequence was 
obtained from ensembl. Oligonucleotide sequences were 
synthesized by Eurofins MWG Operon. 

Primers were designed to include exons and a minimum 
of 50 base pairs (bp) of introns for desired sequence to be 
completely amplified. 

Polymerase chain reaction (PCR) optimization 

The best conditions for PCR were determined using gradient 
PCR. Primers were hydrated with a given volume of dH 2 0 (Table 
2) [15,16] (Acharya et al. 2006), and vortexed for 5-7 secs for 
uniform mixing of primers. Concentration varied between 5 
andio pmol/pl. For a 10 pmol/pl concentration, 90pl of dH 20 
was added to lopl of primers. The ‘master mix’ tubes contained; 
Reddy mix (NEB or Promega): I2.5pl x 14 =i75pl; Control DNA: 
1.5 pi x 14 = 2ipl; dH 20 : 6.5pl x 14 = 91 pi; (Dimethyl sulfoxide) 
DMSO: 2.5 pi x 14 = 35pl; Total 23pl x 14 = 322pl 

The different proportions were multiplied by 14 (see above) 
to account for pipette error. DMSO was added (for F1R1 primers 
only) to prevent formation of spurious secondary structures. 
The Promega (GoTaq® Green Master mix) contains bacterially 
derived Taq polymerase 2 X Green GoTaq® Reaction buffer (pH 
8.5), 400pM dATP, 400pM dGTP, 400pM dCTP, 400pM dTTP 
and 3mM MgCl 2 The New England Biolabs (NEB) (OneTaq quick 
load, 2X MM w/std buffer M0486S) contains dNTPs, MgCl 2 
buffer components and stabilizers. The NEB was used due 
to less smearing compared to Promega. PCR was performed 
under appropriate conditions and according to manufacturer's 
guideline. 


Table 2: Mutations in 5 patients 


Patient ID 

Nucleotide change 

Amino acid change 

rs number 

EXON 

Homozygous or heterozygous 

References 

JH5 

c.l 1 59G>A 

P.E387K 

rs55989760 

(CM980502) 

3 

Homozygous 

[15] 

TC7 

c.230T>C +c.l 103G>A 

P.L77P + 
P.R368H 

CM000136, rs28936414 
(CM000137) 

2,3 

Compound heterozygous 

[16] 

HZ9 

c.l 103G>A 

p.R368H 

rs28936414 

3 

Homozygous 

[16] 

MF10 

c.l 1 59G>A 

P.E387K 

rs55989760 

(CM980502) 

3 

Homozygous 

[15] 

BY17 

c.l 568G>A 

P.R523K 

CM066028 

3 

Heterozygous 

(Acharya et al. 2006) 


Compound heterozygous or homozygous states, their reference SNP (rs) or CM numbers and the exons where they were identified 
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Agarose gel electrophoresis 

1% (w/v) agarose gel concentration was used throughout 
the study. 

Gels were made with lg of agarose powder (molecular 
grade) introduced into 100ml of lx TAE (980ml of dH 20 + 20ml 
of TAE buffer) buffer in a conical flask. After swirling gently, 
agarose-TAE mixture was heated for 2mins until the liquid was 
clear and the powder had dissolved completely. After cooling 
it slightly, lopl SafeView (NBS Biologicals) was added. The 
SafeView dye is safer (less harmful) and stains DNA effectively 
to allow its visualization within the Agarose gel. After swirling 
gently and pouring the gel into the mold, a comb was placed 
in the second groove of the mold. After 20-30 minutes of gel 
setting, the combs were removed and the tray and gel were 
introduced into the electrophoresis tank containing sufficient 1 
x TAE buffer. 5pl appropriately sized DNA ladder (hyperladden 
or loobp #32318) was loaded with pipette into the first well. 
5pl of each amplicon was then loaded into the remaining wells 
carefully. Gels were electrophoresed at 100V voltage, 400Amps 
for 20 - 60 mins. 

After electrophoresis, the gel was visualized to show the 
size (length) of the amplicon. This was carried out by UV- 
transillumination (Fluorchem imager) and an integrated 
camera was used to capture the images (Figure 1). 

PCR amplification 

In order to efficiently amplify the desired target DNA, a 
standard protocol was used in all reactions. Standard 25pl PCR 
reactions contained I2.5pl reddy mix (NEB), 2.5pl DMSO, 7pl 
dH 20 , lpl forward primer, lpl reverse primer (Invitrogen) and 
lpl DNA template. 

After PCR optimization, the right Tm for Primer F1R1 and 
F3R3 were determined as 55°c and 56°c respectively. *Taq DNA 
polymerase extends the target strand at a rate of approximately 
lkb per minute. A negative control (containing lpl of dH 20 
instead of DNA) was also incorporated while setting up PCR 
reactions. 



Figure 1: Gel image of Polymerase chain reaction optimization. 

NEB and Promega stands for the two different 'reddy mixes' used, with forward 
(FI) and reverse (Rl) primers (Exon 2). DNA Hyperladderl on extreme end, NEB on 
left set of bands and Promega on the right. The bands are lined between 1000 and 
1 500bp length. Annealing temperatures 50, 51, 52, 53, 55, 56°c are shown. 
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PCR products purification/cleaning 

PCR plate reactions were purified on Agencourt AMPure 
XP magnetic beads (Beckman Coulter, High Wycombe) using 
a Beckman Coulter Biomek NX robot (Beckman Coulter, High 
Wycombe). Individual PCR amplicons were purified using 
Montage spin columns (Millipore Ltd, Watford) according to 
the manufacturer’s instruction. 

Capillary- electrophoresis based chain-termination 
sanger sequencing 

Taq polymerase enzyme, deoxynucleotide triphosphatases 
(dNTPs) and fluorescently-labelled, chain-terminating 
dideoxynucleotidetriphosphatase (ddNTP) were used in 
standard proportions. Purification of the reaction is carried out 
to eliminate the unincorporated ddNTPs [17]. 

In this study, the 10 pi sequencing reaction (master mix) 
was set up skirted well (PCr) plates with BigDye Terminator V 3.1 
Cycle Sequencing Kit (Applied Biosystems), it contained 0.25 pi 
BigDye V 3 . 1 , 1.875 pi Sequencing buffer, 5 - 375 pl dH 20 ±o. 5 pl, 
i-i.5pl of the AMPure purified amplicons, lpl each of forward 
or reverse primer (in different rows). After mounting on the 
thermocycler (Verti 96 well, AB), the following conditions for 
sequencing were used; 96°c for 2mins, 96°c for losec, 55°c for 
20 secs, 6o°c for 4 mins, 8°c for lomins, 4°c hold for infinity. 

Purification of sequencing reactions 

In order the remove the unincorporated dye terminator, 
Agencourt CleanSEQ paramagnetic bead solution (Beckman 
Coulter) and Beckman Coulter Biomek NX robotics were used. 
Multichannel pipette was used to introduce 5pl of CleanSEQ 
beads into skirted well (PCR tubes containing the sequencing 
product). The purification procedure was performed according 
to manufacturer’s guideline. 

Sanger sequencing analysis 

Two ABI files (sequence traces) for each patient, were 
created and named forward and reverse primer (for example 
JH5 EXON 2F.abi and JH5 EXON 2R.abi respectively). These files, 
and the reference sequence ( CYP1B1) obtained from ensemble 
(ENST00000260630), were uploaded into the GeneScreen 
program. 

The variants (boxes/cells) were edited and the 
chromatogram displayed was observed. Confirmed variants 
occurred in both forward and reverse reads. Strong peaks in the 
chromatogram were only regarded as real sequence changes. 
After analysis, samples with poor results were repeated. 5-10 
nucleotides including variant identified were manually copied 
and matched on the cDNA sequence from ensembl. The number 
(cDNA position) was noted. Reference SNP (rs) number and 
other information such as chromosomal location were obtained 
using ensemble or Exome Variant Server (EVS). 

The staden program helps to efficiently analyse sequences 
and identify variants/ mutations in DNA sequence. The 
Staden Pregap4.i.4bi package was used to align the trace 
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files to the desired reference database. First, the CYP1B1 
(ENST00000260630) gene sequence obtained from UCSC 
genome browser was saved in a text file and uploaded to the 
Staden Pregap4.i.4bi package. Then the Staden Gap package 
helped in analysing the contig (trace sequence) and comparing 
it with a known normal control (cell line DNA). 

Results 

CYP1BJ mutational analysis 

Exons 2 and 3 were (Polymerase chain reaction) PCr- 
amplified and sequenced for each patient as described in the 
Materials and Methods. Sequence files were read with either 
the GeneScreen or Staden programs. One patient DNA (sample 
AKi) failed to give readable sequence (in both exons) despite 
multiple attempts. 13 out of 19 patients had both exons fully 
sequenced (Figure 1). 

In tables 1,3 non-synonymous mutations were identified 
in five probands. In one patient only a single mutation was 
found meaning that the pick-up rate for CYP1B1 mutation being 
pathogenic in this panel of 15 ‘fully sequenced’ patients was 
4/15 (26.67 %). All 3 mutations had been previously reported. 
Four ( 75 %) of these mutations were found in exon 3 while one 
(25%) was identified in exon2. Exon 3 was fully sequenced in 
all 19 patients, while exon 2 was only fully sequenced in 15. 

Out of the five probands that had two mutations, three were 
homozygous for a CYP1B1 mutation and one was compound 
heterozygous. The mutant alleles segregate with the disease in 
an autosomal recessive manner of inheritance [9,18]. Mutation 
c.H59G>A, a nucleotide change at cDNA position 1159 of the 
CYP1B1 gene sequence is the most prevalent homozygous 
mutation found in two patients (Table 2) [15,16] (Acharya et 
al. 2006). 

Single nucleotide polymorphisms (SNP) analyses 

Sequencing of 19 unrelated probands with PCG detected 12 
SNPs which included four intronic variants (two in intron 1 and 
two in intron 2), five variants in exon 2, and 3 variants in exon 
3 (Table 2) [15,16] (Acharya et al. 2006). All SNPS identified in 
introni and protein coding exons 2 and 3 have been previously 
reported in dbSNP, while the changes found in intron 2 are 
novel to the best of our knowledge (Table 2) [15,16] (Acharya 
et al. 2006). Intronic variant c.1043+75 t>g, was the most 
frequent intronic change occurring in 18/19 patients, followed 
by intronic change C.-1-12OT found in 8/19 patients, and c.- 
1-14OT and c.i043+84a>c, each in 1/19 patient. 

The most frequent exon 3variants were c.i294G>C and 
c.i347T>C which were each found in 16/19 patients, followed 
by c.i358A>G which was present in 6/19 patients. 

In exon 2, C.142OG was the most frequent SNP and was 
found in 9/16 patients. This was followed by c.355G>T (7/19) 
and C.777OT, c.729G>C and c.33iG>C (all in 1/19). From the 
result, it could be deduced that patients who had pathogenic 
mutations were homozygous for C.-1-12OT variant. However, 
only three of the mutant patients had this intronic change. 
Mutant patients, who had coding (exonic) SNPs, had them in 


Citation: Ihesiulor GC, Manson F, Ubani UA (2018) A Laboratory Study on the 
DOI: http://doi.org/! 0.17352/2455-1414.000049 


homozygous states except one of the compound heterozygous 
mutant probands (Table 4). 


Table 3: Mutations and SNPs in 5 patients. 


Patient ID 

Variants 

Heterozygous or 
homozygous 

rs or CM 

number 

exon or intron 

TC7 

c.230T>C p.L77P 

Het 

CM000136 

exon 2 

TC7 

c.l 043+75t>g 

Horn 

Novel 

intron 2 

TC7 

c.l 103G>A p.R368H 

Het 

rs28936414 

exon 3 

TC7 

c.l 294G>C p.V432L 

Het 

rsl056836 

exon 3 

TC7 

c.l 347T>C p.D449D 

Het 

rsl056837 

exon 3 

TC7 

c.l 358A>G p.N453S 

Het 

rsl800440 

exon 3 

HZ9 

c.l 043+75t>g 

Horn 

Novel 

intron 2 

HZ9 

c.l 103G>A p.R368H 

Horn 

rs28936414 

exon 3 

MF10 

c.-l-l 2C>T 

Horn 

rs2617266 

intron 1 

MF10 

c. 142C>G p.R48G 

Horn 

rsl 0012 

exon 2 

MF10 

c.l 043+75t>g 

Horn 

Novel 

intron 2 

MF10 

c.l 1 59G>A p.E387K 

Horn 

rs55989760 

exon 3 

MF10 

c.l 294G>C p.V432L 

Horn 

rsl056836 

exon 3 

MF10 

c.l 347T>C p.D449D 

Horn 

rsl056837 

exon 3 

JH5 

c.-l-l 2C>T 

Horn 

Novel 

intron 1 

JH5 

c. 142C>G p.R48G 

Horn 

rsl 0012 

exon 2 

JH5 

c.355G>T p.Al 19S 

Horn 

rsl056827 

exon 2 

JH5 

c.l 043+75t>g 

Horn 

Novel 

exon 2 

JH5 

c.l 1 59G>A p.E387K 

Horn 

rs55989760 

exon 3 

JH5 

c.l 294G>C p.V432L 

Horn 

rsl056836 

exon 3 

JH5 

c.l 347T>C p.D449D 

Horn 

rsl056837 

exon 3 

BY17 

c.l 043+75t>g 

Horn 

Novel 

intron 2 

BY17 

c.l 294G>C p.V432L 

Horn 

rsl056836 

exon 3 

BY17 

c.l 347T>C p.D449D 

Horn 

rsl056837 

exon 3 

BY17 

c.l 568G>A p.R523K 

Het 

CM066028 

exon 3 


Patient ID stands for patient identification, heterozygous (Het) or homozygous 
(Horn), their reference SNP (rs) or CM numbers and the introns or exons where 
they were identified. 


Table 4: The exonic and intronic variants with their rs (reference numbers and 


chromosomal locations. 


Location 

Nucleotide 

change 

Amino acid 
change 

Allele frequency 

(%) 

Reference SNP 

number/ 

INTRON 1 

c.-l-l 2C>T 

NA 

- 

rs2617266 

INTRON 1 

c.-l-l 4C>T 

NA 

- 

rs4987134 

INTRON 2 

c.l 043+75t>g 

NA 

- 

Novel 

INTRON 2 

c.l 043+84a>c 

NA 

- 

Novel 

EXON 2 

c.l 42C>G 

p.R48G 

0=0.336/731 

rsl 0012 

EXON 2 

c.355G>T 

p.Al 19S 

A=0.322/702 

rsl056827 

EXON 2 

c.777C>T 

p.R259= 

NA 

rs72481804 

EXON 2 

c.729G>C 

p.V243= 

G=0.065/142 

rs9341249 

EXON 2 

c.331 G>C 

p.GlllR 

NA 

rs72481804 

EXON 3 

c.l 294G>C 

p.V432L 

0=0.387/843 

rsl056836 

EXON 3 

c.l 347T>C 

p.D449= 

0=0.387/843 

rsl056837 

EXON 3 

c.l 358A>G 

P.N453S 

0=0.099/215 

rsl800440 


*purple coloured variants are novel 
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These four coding SNPs; R48G, A119S, V432L, N453S 
combined as RAVN haplotype, are inherited together as major 
background haplotype for identified mutations [9]. 

Exons are nucleic acid sequences or the expressed region 
of the gene, that codes for a protein. While introns are the 
intragenic regions or intervening sequences that do not code 
for proteins but are removed through (Ribonucleic Acid) RNA 
splicing [19] (Figure 2). 

Variants found in five PCG patients 

Patient JH5: c.H59G>A; P.E387K 

Brief clinical details and family history: Patient JH5 is 
23 years old. Mother may have had congenital glaucoma and 
father has cleft palate (CLP). JH5 presented with symptoms of 
isolated and bilateral glaucoma (PCG). JH5 also has CLP. 

Sequencing analysis shows a previously reported 
homozygous missense change c.H59G>A; P.E387K in exon 
3 of CYP1B1 gene. Five previously reported polymorphisms 
were identified in this patient (c.-i- 12C>T, c.i42C>G p.R48G, 
c. 355G>T P.A119S, c.l294G>C P.V432L, c. 1347T>C P.D449D) 
and one novel variant (c.i043+75t>g) (Figure 3). 

P.E387K is a lysine to glutamine substitution that results 
from a homozygous 1159G>A transition. It occurs in a highly 
conserved region of CYP1B1 exon 3 and affects the conserved 
K helix region of the CYP1B1 molecule. Also, this mutation is a 
founder PCG mutation in the Gypsy subpopulation of Slovakia. 
Cleft palate results from mutations in Transcription factor AP2 
( TFAP2A ) gene. TFAP2A may contribute to anterior segment 
development. 

Patient TC7: c.230T>C P.L77P; c.ll03G>A P.R368H 

Brief clinical details and family history: TC7 is 3 years 
old, born prematurely and presented with bilateral congenital 
glaucoma. 

Two previously reported heterozygous missense changes 
c. 230T>C and c.H03G>A were identified in exons 2 and 3 
respectively. Three coding polymorphisms (P.N453S, P.V432L 
and P.D449D) and one noncoding variant (c.i 043 + 75 t>g) were 
identified (Figures 4,5). 

Patient HZ9: c.ii03G>A; P.R368H 

Brief clinical details and family history: Patient HZ9 is a 3 
years old child of a consanguineous union who has congenital 
glaucoma, dysplastic kidneys and seizures that may be 
syndromic. Also, this patient has done affy array but no result 
was given. 

Sequencing analysis shows previously reported homozygous 
missense change c.H03G>A; P.R368H in CYP1B1 exon 3. One 
novel noncoding variant (c.i043+75t>g) was also identified 
(Figure 6). 

G>A transition at nucleotide 1103 produces an arg368-to- 
his (R368H) mutation found in Saudi Arabian PCG patients. 



Figure 2: Diagrammatic representation of CYP1B1 (Ensembl transcript 
ESNT00000260630) showing the 3 exons, introns and the variants identified. The 
length of exons 1, 2 and 3 are 345, 1044 and 3707 base pairs respectively. The 
coding exons are shaded blue. Mutations are shown in black boxes (above) and 
SNPs below. 
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Figure 3: Staden chromatogram traces of missense mutation c. 1159G>A 
(p.E387K). The top trace shows the mutant sequence and the bottom trace is the 
reference wildtype sequence. Found in patients JH5 and MF10. 
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Figure 4: Staden chromatogram traces of missense mutation c.230T>C (p. L77P) 
in TC7. The top trace shows the mutant sequence and the bottom trace is the 
reference wild type sequence. 


This change showed incomplete penetrance and was reported 
to be absent in 100 Saudi Arabian control chromosomes. Also, 
this mutation has been reported to encode a protein with 
reduced, not abolished enzymatic activity. 

Patient MF10: c.ll59G>A P.E387K 

Brief clinical details and family history: Patient MF10 
is 3 years old and presented with bilateral anterior segment 
dysgenesis syndrome (Peter’s anomaly) manifesting as corneal 
clouding. This condition was present in sister and cousin, and 
the child is of nonconsanguineous parents. 
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Sequencing analysis shows previously reported homozygous 
missense change c.H59G>A P.E387K in CYP1B1 exon 3. Also, 
results show three coding SNPs and 2 noncoding variants 
(Table 4). 

P.E387K mutation is a founder PCG mutation in the Gypsy 
population as previously described in patient JH5. 

Patient BY17: c.1568G>A P.R523K 

Brief clinical details and family history: Patient BY17 is 
9 years old that presented with congenital glaucoma, heart 
murmur delay and Rieger phenotype. This patient has terminal 
deletion of a gene on Chromosome 6 (result not shown). 

Sequencing analysis showed a single heterozygous CYP1B1 
mutation on exon 3. In addition, two coding and one noncoding 
variant were found in this patient (Figure 7). 

Congenital glaucoma in this patient was not due to CYP1B1 
gene mutation. As stated earlier, glaucoma-associated Axen- 
field Rieger syndrome is due to deletion of gene located in 
chromosome 6p25 (forkhead Box Cl gene). This is responsible 
for the glaucoma in this patient. 
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Figure 5: Staden chromatogram traces showing missense mutation c,1103G>A 
(p.R368H) in a heterozygous state. The top trace shows the mutant sequence and 
the bottom trace is the reference wildtype sequence. Found in TC7. 
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Figure 6: Staden chromatogram traces of missense mutation c,1103G>A 
(p.R368H) in a homozygous state. The top trace shows the mutant sequence and 
the bottom trace is the reference wildtype sequence. Found in patient HZ9. 



Figure 7: GeneScreen chromatogram traces showing missense mutation 
c.l 568G>A (p.R523T). The top trace shows the forward sequence which has the 
mutation (c.l 568G>A) in BY17, and the bottom trace is the reverse sequence. 


Table 4: The exonic and intronic variants with their rs (reference numbers and 
chromosomal locations. 


Location 

Nucleotide 

change 

Amino acid 
change 

Allele frequency 

(%) 

Reference SNP 

number/ 

INTRON 1 

c.-l-l 2C>T 

NA 

- 

rs2617266 

INTRON 1 

c.-l-14C>T 

NA 

- 

rs4987134 

INTRON 2 

c.l 043+75t>g 

NA 

- 

Novel 

INTRON 2 

c.l 043+84a>c 

NA 

_-_ 

Novel 

EXON 2 

c.l 42C>G 

p.R48G 

0=0.336/731 

rsl 0012 

EXON 2 

c.355G>T 

P.A119S 

A=0.322/702 

rsl 056827 

EXON 2 

c.777C>T 

p.R259= 

NA 

rs72481804 

EXON 2 

c.729G>C 

p.V243= 

G=0.065/142 

rs9341249 

EXON 2 

c.331 G>C 

p.GlllR 

NA 

rs72481804 

EXON 3 

c.l 294G>C 

p.V432L 

C=0.387/843 

rsl 056836 

EXON 3 

c.l 347T>C 

p.D449= 

C=0.387/843 

rsl 056837 

EXON 3 

c.l 358A>G 

P.N453S 

C=0.099/215 

rsl800440 


*purple coloured variants are novel 


Discussion 

CYP1B1 gene mutations are the major molecular cause of 
PCG [20]. PCG is a genetically heterogeneous phenotype [13]. 
In this research, the genetic analysis of the coding exons 2 & 
3 in CYP1B1 gene (Ensembl transcript ESNT00000260630) of 
15 unrelated PCG cases is carried out. It shows that about 4 
(26.67%) out of 13 probands had at least one PCG-causing 
mutation previously identified in scientific research. In the four 
patients having mutations, three had homozygous mutations 
while one showed compound heterozygosity. 

In addition, one patient (BY17) had a single heterozygous 
mutation in CYP1B1 and a terminal deletion in forkhead 
transcription factor gene ( FKHL 7). The deletion in the short 
arm of chromosome 6 (6p25) is responsible for Axenfeld- 
Rieger Syndrome (ARS) and glaucoma [21]. Patient (TC7) was 
born prematurely. And it is known that anterior segment 
anomalies can be present in premature babies [22,23]. 

The pick-up rate (26.67%) may be because sequencing 
was limited to CYP1B1 coding exons leaving the promoter or 
non-coding regions [8]. Other reasons could be that the PCG 
phenotype may be caused by mutation in other PCG-associated 
loci (such as GLC3B, GLC3C or GLC3D ( LTBP2 ) [8]. Pathogenic 
mutation in the LTBP2 gene associated with PCG has been 
reported in Pakistani, European Gypsy and Iranian probands 
[24,25]. Moreover, it has been suggested that a pathogenic 
mutation in MYOC gene may cause PCG with a CYP1B1 mutation 
via a digenic mode [26]. However, this finding is not established. 
Mutation in other unknown genes may be responsible for PCG 
in the CYP1B1 negative patients [8]. 

The rest fourteen patients were CYP1B1 negative. BY17 
has ARS, so can be discounted. Although exons 2 in four 
probands (JH5, MF10, CS15 and JL21) was not fully sequenced, 
homozygous pathogenic mutations were found in two (JH5 and 
MF 10) (Table 2) [15,16] (Acharya et al. 2006). 
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From our results, PCG showed allelic heterogeneity in the 
patients (JH5 and MF10) associated with homozygosity for 
2 different CYP1B1 mutations and on patient (TC7) showing 
compound heterozygosity in two distinct CYP1B1 mutations. 
Allelic heterogeneity explains the molecular contribution for a 
uniform clinical manifestation of probands, and homozygous 
and compound heterozygous explains the autosomal recessive 
nature of PCG [15,27]. 

P.R368H (found in TC7 and HZ9) is a hypomorphic mutation 
located in the conserved core structure (CCS), J helix region 
of the CYP1B1 protein [16]. This causes a decline in CYP1B1 
enzymatic activity (~ 20% of wild type protein) or decrease in 
protein stability [27]. In a study by Li et al. [20], P.R368H was 
mostly found in white Europeans [20]. 

Mutation P.E387K (found in JH5 and MF10) is most 
common among Slovakian Gypsies, less in Caucasians (4.90%), 
present in Amish population and absent in Asians and middle 
easterners [18,20]. For effective genetic screening in PCG 
patients, it is necessary to find out the common and founder 
mutations in a given population [21]. 

E387K is the founder mutation that accounts for 79.63% 
of CYP1B1 gene mutations [15]. The high rate of consanguinity 
(especially cousin to cousin marriages) and high coefficient 
of inbreeding in the Middle Eastern and Gypsy populations 
explains the high occurrence of the E387K mutation [20]. 

The amino acid position (GIU387) in which P.E387K 
mutation changes glutamic acid to lysine, is a highly conserved 
position in all documented species and P450 enzymes [28]. 
Also it’s a core element located in helix K, a region suspected to 
be essential for proper folding and active haem binding of the 
CYP1B1 enzyme [19,29]. 

The mutation found in patient BY17 (c. 1568G>A) was 
previously identified in Israeli Bedouin kindred and it was 
stated that it obliterates the Ddel restriction site and disrupts 
the CYP enzyme active site found in the C-terminal region of 
the CYP1B1 enzyme [27]. Bar-yoseph et al. [27], reports this 
variant that changes arginine to lysine in amino acid position 
523 (P.R523K), was not found in 100 healthy individuals [27]. 

In our study, fourteen patients had no disease-associated 
mutations in the fully sequenced regions of the CYP1B1 gene. 
This supports the report that some unidentified molecular 
etiology is behind the PCG in some patients [27]. 

Residue change from leucine (L) to Proline (P) at position 
77 (P.L77P), results from c.230T>C missense mutation. This 
previously reported mutation (found in patient TC7) occurs in 
conserved position. It has been reported to be associated with 
PCG in Saudi families [16]. 

Information from this study shows that early genetic testing 
is pertinent to determine the carrier status of individuals and 
their phenotype [20]. Late presentation of PCG is associated 
with profound visual impairment in children [21]. 

Having identified CYP1B1 negative patients, future studies 
can be undertaken to determine other genetic causes of PCG by 
whole exome sequencing of these subjects. 
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Furthermore, the absence of CYP1B1 mutations in some PCG 
patients supports that another unidentified gene is mutated. 
The residual level of CYP1B1 activity is modulated by the 
presence of modifier genes [9]. 

Conclusions 

The number of CYP1B1 negative (11/15) probands was higher 
than the CYP1B1 positive (4/15) patients. Some mutations may 
have been missed due to the genetic screening strategy applied. 
This is the first drawback of this study. Screening was not 
extended up to the regulatory sequences in the upstream or 
downstream regions in the introns although bit of introns were 
sequenced due to primer design [9]. 

The noncoding regions (including exon 1) may have the 
disease-related variants whereas this study screened only 
protein coding exons 2 and 3 of the gene [8]. For accurate 
genetic diagnosis of PCG, all three exons need to be sequenced 
[ 9 ]. 

Secondly, the small sample size (19) prevents conclusions 
to be confidently made; confirmation of result is therefore 
needed with larger number of subjects [20]. 

In conclusion, the analysis of CYP1B1 genotypes of 19 
unrelated patients is reported in this study. A relatively low 
proportion of our subjects (~ 27%) tested positive for PCG- 
related CYP1B1 mutations demonstrate the need to identify 
other PCG-causing genes. 

Exome sequencing of CYP1B1 negative patients to detect 
new genes mutated in PCG. 

Co-segregation analysis of other unaffected and affected 
family members may be carried in the future. 
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